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. Aims. High-resolution spectroscopic observations (UV HST/STIS and optical) are used to characterize the physical 

state and velocity structure of the multiphase interstellar medium seen towards the nearby (170 pc) star HD 102065. 
The star is located behind the tail of a cometary-shaped, infrared cirrus-cloud, in the area of interaction between the 
, Sco-Cen OB association and the Local Bubble. 

Methods. We analyze interstellar components present along the line of sight by fitting multiple transitions from a group 
of species all at once. We identify four groups of species: (1) molecules (CO, CH, CH + ), (2) atoms (C I, S I, Fe i) with 
ionization potential lower than H I, (3) neutral and low-ionized states of atoms (Mg I, Mg II, Mn II, P II, Ni II, C II, 
N I and O i) with ionization potential larger than H I and (4) highly- ionized atoms (Si ill, C IV, Si iv). The absorption 
spectra are complemented by H I, CO and C II emission-line spectra, H2 column-densities derived from FUSE spectra, 
?— 1 ■ and IRAS images. 

' Results. Gas components of a wide range of temperatures and ionization states are detected along the line of sight. Most 

of the hydrogen column-density is in cold, diffuse, molecular gas at low LSR velocity. This gas is mixed with traces of 
warmer molecular gas traced by H2 in the J > 2 levels, in which the observed CH + must be formed. 
We also identify three distinct components of warm gas at negative velocities down to — 20kms _1 . The temperature 
and gas excitation are shown to increase with increasing velocity shift from the bulk of the gas. 
Hot gas at temperatures of several 10 5 K is detected in the most negative velocity component in the highly-ionized 
specie. This hot gas is also detected in very strong lines of less-ionized species (Mg 11, Si 11* and C 11*) for which the 
bulk of the gas is cooler. 

Conclusions. We relate the observational results to evidence for dynamical impact of the Sco-Cen stellar association on 
the nearby interstellar medium. We propose a scenario where the infrared cirrus cloud has been hit a few 10 5 yr ago by 

, a supernova blast wave originating from the Lower Centaurus Crux group of the Sco-Cen association. 

■ The observations provide detailed information on the interplay between ISM phases in relation with the origin of the 

' Local and Loop I bubbles. 

Key words. ISM:structure, ISM:clouds, ISM:kinematics and dynamics, ISM:individual objects:Chamaeleon clouds, 
Ultraviolet:ISM, Stars:individual:HD102065 



■ 1. Introduction We hereafter refer to it as the Blue Cloud in reference to 

its blue infrared colors. While the main Chamaeleon clouds 
Astronomical observations provide multiple perspectives on have a com p ar able extent in CO surveys as in IRAS im- 
the interstellar matter, gas and dust, that are rarely gath- ages , only the center of the Blue Cloud head is detected in 
ered together to study the relation between interstellar- co dBoulanger et al. 19981 : iMiluno et al. 200lh . Boulanger 
medium phases. It is the ambition of this paper to present et aL estimate that only 10% of the cloud mass is seen 
such a study on gas and dust observed in the direction of the in co emission. On larger angular scales, HD 102065 lies 
nearby (170 pc), lightly reddened (E(B-V) = 0.17), B9IV within the sky area where various evideil ces for an inter- 
star HD 102065 m Chamaeleon (Table 1). In IRAS images, action between the Scorpius-Centaurus (Sco-Cen) OB as- 
the infrared cirrus (Dcld 300.2-16.9) seen in the foreground soc i at i on and the Local Bubble have been gathered from 
of HD 102065 has a prominent cometary shape and an un- H t observ ations (de Geus 1992), the X-ray ROSAT survey 
usually blue IRAS colors (high 12 and 25 ^m to 100 ( Egger and Aschenbach 1995) and optical absorption-line 
brightness ratio), indicative of a large abundance of small spectroscopy (Corradi et al 2004) 

stochastically-heated dust particles (Boulanger et al. 1990) . HD 102065 is located at the same Galactic longitude and 

slightly lower Galactic latitude than the Lower Centaurus 



Correspondence to: cecile.gry@oamp.fr Crux (LCC) group of the Sco-Cen association. Hipparcos 



positions, proper motions, and parallaxes have been used to 
show that the LCC and the Upper Centaurus Lupus (UCL) 
group of stars, presently located within the Loop I bubble 
at a distance of ~ 120 pc, were located closer to the Sun 
5 — 7Myr ago. In two independent works (Maiz-Apellaniz 
2001 and Bcrghofer and Breitschwerdt 2002), this informa- 
tion was used to relate the origin of the Local Bubble to the 
~ 20 supernovae explosions associated with the LCC and 
UCL groups that are estimated to have occurred in the 
last 10 Myr. The structure and physical state of the inter- 
stellar medium seen in direction of the Sco-Cen association 
should reflect the action of the supernovae explosions pow- 
ering and progressively shredding, and spreading away the 
stars' parent cloud. This scenario has been introduced into 
a numerical simulation based on the supernovae-drivcn ISM 
model of de Avillez (2000) by Breitschwerdt and de Avillcz 
(2006). In the reproduction of the local interstellar medium 
generated by this simulation (Figure 1 in Breitschwerdt and 
de Avillez), HD 102065 is located within the Loop I bub- 
ble and the line of sight to the star crosses the shell that 
provides separation from the Local Bubble. 

Observations characterizing the gas and the dust along 
the line of sight to HD102065 have been presented in ear- 
lier publications (Boulangcr, Prevot & Gry, 1994 ; Gry 
et al, 1998). Unlike most early- type stars observed in UV 
spectroscopy, the position of HD 102065 corresponds to no 
enhancement in IRAS emission, implying that the mat- 
ter responsible for ab sorption is not heated by the star 
(|Boulanger et al. 19941 ). The extinction curve produced us- 
ing IUE low-resolution spectroscop y, shows a strong 220 nm 
bump and a weak far-UV rise (|Boulanger et al. 19941 ). 
Moderate-resolution spectroscopy (A/AA = 20 000) ac- 
quired by FUSE, detected molecular hydrogen absorption 
lines. These observations showed that H2 at a temperature 
of ~ 60 K is the dominant state of hydrogen along the line 
of sight to HD 102065. The presence of a smaller amount 
of war mer H2 (a few 10 K) was inferred from the H2 exci- 
tation dGrv et al. 20021) . HD 102065 was also observed us- 
ing the Goddard High Resolution Spectrometer (GHRS) 
on the Hubble Space Telescope (HST). These observations 
revealed the presence of high, negative-velocity gas with 
unusually high Si ii excitation that was seen as the signa- 
ture of the dissip ation of a large amount of kinetic energy 
dGrv et al. 19981) . 

A number of questions remain unanswered. Why is the 
small dust abundance enhanced in the Blue Cloud, which, 
from the FUSE and optical spectroscopy, appears to be 
a typical diffuse molecular cloud? How was the negative- 
velocity gas accelerated? Are the low and negative-velocity 
gas physically connected? More generally, how can the data 
be interpreted within the present understanding of the in- 
teraction between the Loop I and the Local bubbles? 

In this paper, we present high-resolution (R = 100 000) 
UV spectr,a obtained using the HST Imaging Spectrograph 
(STIS) complemented by optical absorption spectra, and 
CO and H 1 emission spectroscopy (Sections[2]and 3). These 
data improve previous observations in spectral resolution 
and spectral coverage. The multi-wavelength spectra are 
used to describe the physical state and velocity structure of 
the multiple gas components observed along the line of sight 
(Section H}. In Sectional we discuss the cool, low- velocity 
gas, and in Section we relate the observational results 
to evidence for dynamical interaction between the Sco-Cen 
stars and the nearby interstellar medium, and discuss the 
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Table 1. HD 102065 "fact shee t" . A v is the visible ex tinc- 
tion, E(B-V) is the col or index dBoulanger et al. 1994f ). H 2 
column densities from (jGrv et al. 20021 )! The total column 
density (Njj = N(H) + 2 N(H 2 )) has been derived from 
E(B-V). 



possible interplay between the multiple components. The 
main conclusions of this paper are gathered in Section [7] 

2. Ultraviolet absorption lines: STIS data. 

2.1. Observations 

We acquired STIS UV spectra of HD 102065 with the high- 
resolution MAMA Echelle gratings E140H and E230H. We 
used five different grating settings covering the following 
wavelength regions: 1144 - 1324 A, 1322 - 1512 A, and 1494 
- 1674 A (E140H with aperture 0.2 x 0.09), and 2376 - 2640 
A and 2620 - 2864 A (E230H with aperture 0.1 x 0.09). 

We performed spectral extraction using the calibration 
software CALSTIS version 2.14c. This version provided ro- 
bust calibration of the zero-flux level, an important param- 
eter when deriving column densities from absorption lines: 
in all strongly saturated lines, the mean zero-level is 4 to 5 
times below the RMS noise. We therefore assume that the 
zero-flux level is valid for all spectra, and we do not consider 
it a free parameter when fitting the absorption profiles. 

For the wavelength region 1144 - 1324 A, observations 
were completed in five sub-exposures. These five measure- 
ments do not have the same wavelength grid and for each 
wavelength step the five wavelength values are distributed 
over the entire resolution element. In order to not degrade 
the resolution when co- adding the five sub-exposures, we 
have defined a mean wavelength grid and have interpolated 
each sub-exposure spectrum over this grid in common. The 
resulting smoothing factor with this method is no more 
than 15%. 

The resolution achieved with the echelle gratings is 
about 2.6kms _1 . The line spread functions (LSF) are tab- 
ulated in the STIS Instrument Handbook. We have per- 
formed a double-Gaussian fit to the tabulated LSF and use 
the resulting FWHM values in the profile-fitting software. 
The results of the double-Gaussian fit to the LSF were a 
narrow component with FWHM ~ 1.2 pixels (6.3 mA) for 
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Table 2. Wavelength and /-values for the observed absorption lines. Most of those values are from Morton (2000 and 
2001) and Verner et al. (1994). Oscillator strengths values for S I are taken from Beideck et al. (1994). Values for O I 
1355.5977 A and Ni II 1317.2170 A are from Welty et al. (1999). CO oscillator strengths have been kindly provided by 
Dr. Eidelsberg. 



E140H at 1200 A, 1.1 pixels (7.2 mA) for E140H at 1500 A, 
and 1.7 pixels (17.9mA) for E230H at 2400 A, and a broad 
component with FWHM ~ 6 pixels in all cases. 

According to the STIS Instrument Handbook, the abso- 
lute wavelength accuracy across exposures is 0.2-0.5 pixels. 
We verified that these measurements were applicable to our 
own observations by comparing the centers of spectral lines 
corresponding to the same element, occurring in different 
exposures at different wavelengths: shifts between differ- 
ent wavelength ranges were never larger than 0.5 pixel, or 
0.6 km s -1 . 



The achieved signal-to-noise ratio (hereafter SNR) in 
the continuum, ranges from 23 to 41 depending on the grat- 
ing setting. At the positions of most interstellar lines the 
SNR ranges from 20 to 40, except for the lines Si II*, C n 
and C II*, where it drops to 6 because they are located at 
the bottom of strong stellar lines. 
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Figure 1. Molecular lines. From top to bottom: 12 CO J=2- 
1 emission line (Section 3.3) in the millimeter range with the 
SEST, two examples of the CO ultraviolet absorption bands 
(A-X (2,0) and A-X (0,0)) observed with STIS, and the CH 
and CH + visible absorption lines (Section 3.1). Histogram- 
style curves represent the observations normalized to the 
continuum, solid lines are the best fits, and dashed lines, the 
profiles before convolution with the instrument line spread 
function, i.e. the intrinsic interstellar profile. The letters A, 
B and C indicate velocity components that will be intro- 
duced in the spectral analysis. 



2.2. A guided tour of the STIS absorption spectra 

Here we provide an overview of the interstellar, absorption 
lines present in the STIS spectra and used in the analysis. 
They are listed in Table [2] 

We identify four groups of species, classified according 
to the characteristics of their line profiles, which reflects 
the distribution and the nature of the gas traced by the 
given species. Examples of the line profiles are presented in 
Figures [T] to |U The figures also show the velocity compo- 
nents and fits to the profiles, which are the results of the 
analysis which will be described in Section [U after presen- 
tation of the complementary data in Section [31 

(f) Molecules: see Figure [TJ Several CO electronic A- 
X bands are present in the HST/ST1S spectra. We use the 
four following bands where CO is detected in several lines 
from rotational levels J"=0, 1 and 2 in the HST/STIS ob- 




1425.00 14:35.05 1425.10 1425.15 1425.20 1279.85 1279.90 1270.95 1280.00 1280.05 



SI 1295. Go A 



1295.60 12S5 65 1295.70 1295.7 




CI* 1280.59 A 



1280.55 1280.60 12S0 65 1230.70 1280.75 



6xlU" 13 Ej- 



. Fcl 2484.02 A 




-ts . Cf 1329.6 A 
CI" 1329.6 A 



1329.55 1329. GO 1329 65 1:329.70 1329.75 




CI 1276. 4R A 



1276.45 1276.50 1276. 

Wavelength A 



1276.60 1277.50 1277.55 1277.60 1277.65 1277.70 

Wavelength A 



Figure 2. Sample lines from neutral species. Histogram- 
style curves represent the observations (in erg cm -2 s~ 1_1 ), 
solid lines are the best fits, and dashed lines, the profiles 
before convolution with the instrument line spread func- 
tion, i.e. the intrinsic interstellar profile. An asymmetry is 
evident on most profiles which arc therefore best fitted with 
two components, called A and B. 



servations : A-X (0,0) at 1544.54 A, A-X (1,0) at 1509.9 A, 
A-X (2,0) at 1477.57 A and A-X (3,0) at 1447.4 A. Several 
other A-X bands are present in the spectra but they are 
both weaker and noisier (due to their location at shorter 
wavelengths hence in a fainter part of the spectrum). The 
bands are resolved into R and Q lines from the three ro- 
tational levels J" = 0, 1 and 2, as illustrated in Figure [1] 
where two bands are shown as examples. 

(2) Neutral species with ionization potential lower 
than that of H I: see Figure [5] 

Atomic carbon lines have been measured in the three fine- 
structure levels of the 3 P fundamental level: J = (here- 
after C i), J = 1 (hereafter C I*), and J = 2 (hereafter 
C I**) . 36 carbon lines have been detected, among which we 
have selected 20 for use in the fits, eliminating lines which 
are either heavily blended with another line, or lack reliable 
f-values. A representative subset of those lines is shown in 
Figure [21 with lines from the other neutral species, Fe and 
S. These neutral lines are narrow with a width of only a few 
resolution elements (i.e. a few kms -1 ), and clearly asym- 
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Figure 3. Sample of lines from ionized species and neutral species with ionization potential close to that of H i. Same 
line style as Figure [5] Most fits have been performed with the same line-of-sight model, made of 4 components with 
same velocity and same temperature for all species. There are only three exceptions: (i) An additional component of 
high-temperature gas is included at the velocity of component 4 to fit the extended blue wing of the strong Mg n and 
C II* lines, (ii) In the O I** profile an additional component is considered, blueward of component 4. (hi) The S II profiles 
include an extra absorption on the red side of component 1. 
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3.1. Optical absorption lines 

Optical spectra of the CH and CH + lines about 4300 A 
(see Table E]), the C 2 A 1 ^ - X 1 E+(2,0) and CN A 2 n - 
X 2 £+(2,0) molecular bands about 8760 and 7905 A, 
were obtained in April 1990 using the Coude Echcllc 
Spectrometer at the 3.6 m of ESO. 
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Table 3. Molecular absorption lines detected in the optical 
spectra. Listed parameters are from Gredcl ct al. (1993). 
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Figure 4. Highly-ionized species observed with STIS. As 
an indication, consistent absorption line profiles for gas at 
a temperature of 900 000K and at the velocity of component 
4 are superimposed. 



metric indicating the presence of two velocity components 
referred to as A and B. 

(3) Low ionization species and neutrals with ion- 
ization potential higher than that of H I: see FigureEl 
In those species the absorption extends to negative veloci- 
ties as low as — SOkms" 1 for some of the lines. The broad 
absorption is resolved into four distinct components. They 
are particularly wcll-rcsolved in the lines of the heavy cle- 
ment Fe II. These components are numbered 1 to 4 in order 
of decreasing velocity. 

(4) High ionization species: see Figure |H 

High ionization species Si in, Si iv, C iv are detected in this 
line of sight. Their profile is broad, and centered at negative 
velocity close to the most negative velocity component seen 
in the other ions, Component 4. 



3. Complementary observations 

We present additional data sets that complement our STIS 
observations. 



The spectrograph was connected to the 3.6 m beam by a 
fiber link. The light coming from the fiber was divided in 10 
independent slices. Each of the slices produced a 20-pixel 
wide spectrum on the CCD. To reduce the read-out noise 
CCD pixels were read in bins of 10 pixels, perpendicular 
to the dispersion direction. The optics of the spectrome- 
ter were oriented so that the spectral resolution was not 
reduced (i.e. the dispersion direction was perpendicular to 
the binning direction). The resolving power was 110 000 
corresponding to a velocity resolution of 2.7 km s -1 . The 
sampling was 2.2 pixels per resolution element, and the 
spectra were 1024 pixels in length. Wavelength calibration 
was performed using an observation of a Th-Ar cathode 
lamp. This calibration was repeated for each set-up of the 
spectrometer. The identification and line-center measure- 
ments of a set of around 20 lines in these spectra allowed 
us to derive the dispersion relation with an accuracy better 
than 10 mA. Flat-field images were acquired at the be- 
ginning and end of each set of observations, for the entire 
wavelength range. 

The images were corrected for dark current, and then 
flat-field corrected. Spectra were extracted from the images 
by averaging 20 lines over which the flux is spread. The 
spectra were normalized by dividing the original spectra 
by a median average showing no spectral lines from inter- 
stellar gas or the Earth atmosphere. Atmospheric lines were 
identified using observations of HD 106911, a bright star in 
Chamaeleon with almost no reddening. 
Sections of the normalized spectra that include the CH 
and CH + lines, are presented in Figure [TJ The C2 and 
CN spectra show no detections with 3cr upper limits of 
2.7 and 1.0 mA, respectively. Gaussian fits to the CH and 
CH + lines have been used to derive the line parameters, 
and thereby the column density N. Because the transitions 
are optically-thin, N is related to the measured equivalent 
width, W x , by the equation: N = 1.13 10 20 W\/X 2 / with 
W\ and A in A, AT in cm and where / is the line oscillator 
strength given in Table [3] Results are presented in Table El 
The optical absorption lines are also fitted, together with 
the UV absorption lines, in the general spectral analysis de- 
scribed in Section 14.2.21 This, in particular, allows the CH 
and CH + velocity distribution to be related to the velocity 
structure seen in other species. 



3.2. Hi emission 



A sensitive 21-cm H I southern sky survey has recently been 
published (Bajaja et al. 2005). Its angular resolution is half 
a degree and the velocity resolution 1.3 kms -1 . In Figure[5j 
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Figure 5. H I emission line in the direction of HD 102065 
derived from the IAR Southern Sky survey. The thick line 
represents the difference with respect to a reference position 
outside of the Blue cloud. 

we present the survey spectrum for the position closest to 
HD 102065 (thin line). The thick line is the difference be- 
tween this spectrum and an OFF cloud spectrum obtained 
by averaging 3 spectra outside the Blue Cloud (Dcld 300.2- 
16.9) at 1 = 297° to 297.5° and b = -19.5° to -19°, corre- 
sponding to a local minimum in the 100/im IRAS map. This 
subtraction is supposed to remove the background emis- 
sion, but may also remove foreground emission that is not 
spatially correlated with the Blue Cloud. Within unknown 
variations of the background emission between the ON and 
OFF cloud positions, this difference spectrum represents 
the H I spectrum of the Blue Cloud. In Tabled we list the 
results of a Gaussian decomposition of this spectrum. The 
column densities have been calculated with the assumption 
that the H I emission line is optically-thin. 

In Figure [H we show the spatial distribution of the H I 
emission in two-channel maps (Vlsr — — 1 an( i +5 km s~ 1 ) 
overlaid over the IRAS 12/im image (Mivillc-Deschenes and 
Lagachc, 2005), and in Figurc[7]a position-velocity diagram 
across the Blue cloud at the Galactic latitude of HD 102065. 
These two figures show that the Blue Cloud, the bright 
IRAS 12 fim emission on the map, is associated with gas 
at low velocities. The total brightness spectrum (thin line 
in Figure [5|) shows the presence of a broad, negative ve- 
locity component, extending down to —40 kms . This gas 
is not spatially correlated with the infrared emission from 
the Blue Cloud. It is however observed in UV-absorption 
spectra, and must therefore be located in front of the star. 

3.3. CO emission spectra 

Long integrations have been performed with the 15m SEST 
telescope on the 12 CO J=2-l and J=l-0 emission lines. 
The J =1-0 and 2-1 lines were observed in 1990 and 2001, 
with system temperatures of 700 K and 200 K respectively, 




Figure 6. IRAS 12/im image (IRIS processing, Miville- 
Dcschenes and Lagache, 2005) with contours of H I emis- 
sion at Vlsr = — lkms -1 (white dashed) and 5kms~ 1 
(black solid) taken from the IAR survey (30' resolution and 
1.3 kms -1 width velocity channels, Bajaja et al. 2005). The 
position of HD 102065 and the T Tauri star T Cha arc 
marked. Contours are labeled with the brightness tempera- 
ture. The 30 K contour is thicker. The cloud tail is roughly 
at a constant Galactic longitude. The Galactic coordinates 
of HD 102065 are I = 300.0° and b = -18.0°. The elec- 
tronic version is in color. 
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Figure 7. H I position- velocity diagram across the Blue 
Cloud at the Galactic latitude of HD 102065 (b=-18°). This 
figure shows that the two velocity components at respec- 
tively 0.1 and 4.2 kms -1 each extends on one side of the 
star position (1=300°) where they overlap. The emission at 
the OFF position (see text) has been subtracted to high- 
light the gas associated with the cloud. 



and a velocity resolution of 0.1 kms 1 . The rms noise level 
of the 12 CO J =2-1 line displayed in Figure [T] is suffi- 



Sp6C16S 


Vlsr I^Kms ) 


lb (^J 


r VV X11V1 I Kill o 1 


K firm q~M 
D I Kill a 1 


IN I LII1 1 


Hi 21 cm (difference) 


14.5 ± O.O 


a a _i_ a 1 k 
u.y ± U.lO 


f.O ± 1.4 


A A _i_ a o 

4.4 zb U.o 


i o _i_ a o 1 
l.o zh U.z 1U 




4.z ±0.1 


ZU.4 zt 1.1 


o o i a o 
O.Z ± U.Z 


i a _i_ a i 
i.y it U.l 


l.o ± U.U/ 1U 




0.1 ±0.1 


13.4 zb U.7 


4.3 ± 0.3 


2.6 =b 0.2 


i i _i_ a c\a i a20 

1.1 ± U.Uo 1U 




— 1 / .0 It U.D 


l.U zt u.u / 


1 C _L 1 A 
10 zt 1.4 


Q A -4- H Si 

y.4 zt u.o 


o.u zt u.z iu 




no n ri 
Z.O ± U.Z 


7 a i a - ^ 
/ .U 1U 


O.O zt U.o 


o i i a o 
z . 1 zb U . z 


o.y iu ( j = ij 


C0(2-l) 


2 75 + 1 7 


3.6 10~ 2 


3.0 ± 0.3 


1.8 ± 0.2 


1 8 10 13 ( 7 — 2) 






/ (erg s cm sr ) 






MfP TT*^ 

IN 11 J 






Z.O 1U 






i o i n ib 

± .z ±u 



Table 4. Results from emission line measurements. H I measurements are performed on the difference spectrum (see 
text and Figured]). The H I component at 14.5 km s -1 is not detected in absorption and is probably due to background 
gas situated behind the star HD 102065. 



ciently low (~ 10 mK) to reveal a broad and weak compo- 
nent (Aw = 2.9 km s -1 , T* A = 36 771-ftT), and a marginally- 
detected narrower component (Aw = 0.5 km s -1 , T* A = 
l5mK). The antenna temperature has not been corrected 
for beam coupling to the sky, because the source is much 
more extended than the tclcscopc-beam size at this fre- 
quency (~ 22" or 0.015 pc at the distance of the cloud). 
This is not a critical issue because the beam efficiency was 
~ 0.85 at 230 GHz. Similarly, the rms noise level of the 
J =1-0 line is 18 mK for a velocity resolution of 0.2 km s" 1 . 
Given the weakness of the 12 CO lines, the assumption of 
optical thinness is justified and the column densities of 
molecules in the J = 2 and 1 levels are inferred from the line 
integrated intensities W(CO) = 0.1 and 0.008 K kms -1 of 
the broad and narrow components of the J = 2 — 1 line, and 
W{CO) = 0.2 A'kms- 1 for the J=l-0 line. The respective 
column densities are given in Table |U 

3.4. C + observations 

The C + fine-structure transition at 158 /im has been mea- 
sured towards HD 102065 with the Long Wavelength 
Spectrometer (LWS) of the Infrared Space Observatory 
(ISO), with a spectral resolution of R = 300, much lower 
than all other data presented in this paper. The data 
have been processed following the LWS handbook (Gry 
et al 2003). The resulting line intensity I(C + A158^m) is 
2.8 ± 1.0 10 -6 ergs -1 cm -2 sr" 1 , from which we estimate 
the C II* column density N(C II*) = 1.2 ± 0.4 10 15 cm" 2 . 
This is a measurement of the mean column density over the 
80" LWS beam. 

4. Analysis of the line of sight - Description of the 
interstellar components 

4.1. The tool: Multi-lines fitting 

To derive the characteristics of the interstellar components, 
i.e. column density N, velocity v and broadening parame- 
ter b, we fit the observed line-profiles with theoretical ones, 
which are the results of the convolution of a Voigt profile 
with the instrumental function (LSF). We use a profile- 
fitting method, instead of a traditional measurement of 
equivalent width, because the velocity structure is complex 
as can be seen in Figures [U to [31 Because several com- 
ponents are blended together in most lines, profile fitting 
is required to derive the characteristics of the individual 
components. The data fitting has been performed with the 



use of the profile fitting software "Owens. f" developed by 
Martin Lemoine. 

The software allows several lines of the same element, as 
well as lines from different elements, to be fitted together. 
It also allows several velocity components to be fitted si- 
multaneously, and the components to be included in the 
absorption profiles for each element to be specified. The 
software therefore derives a global and consistent solution 
for a group of species. Multiple iterations are necessary to 
find which species could be fitted in a consistent way with 
the same velocity components. A common component for 
different species obviously implies a common absorption ve- 
locity but also common physical conditions, implying con- 
sistent broadening parameters. 

The fitting software breaks the line-broadening param- 
eter (b-value) into thermal broadening, which depends on 
the element mass, and non-thermal broadening (turbu- 
lence), which is the same for all elements in a component. 
Therefore, fitting lines from elements of different masses 
simultaneously, in principle enables the simultaneous mea- 
surement of temperature, turbulent velocity, the velocity 
itself, and column densities. In practice interstellar absorp- 
tion lines are usually found to be of two types. The first 
type displays unsaturated lines from heavy elements that 
have similar b- values independently of the element mass: in 
this case, the broadening is primarily non-thermal, imply- 
ing that T is low and can be estimated only if the analysis 
includes an element whose mass is low enough to introduce 
a significantly different b-value. Here we will consider H I 
emission lines. In the case of the second type of interstel- 
lar absorption line, the b-values vary with element mass 
due to a significant thermal effect, implying that T is high 
and can be estimated. In practice the interstellar absorp- 
tion lines are broader than the resolution element (for T 
higher than a few thousand K), and strong lines are occa- 
sionally broadened considerably by saturation, leading to 
the blending of components and limiting the precision with 
which broadening parameters can be derived. 

The quality of the fits, and the error bars on column 
densities, have been computed using the A% 2 method de- 
scribed, for example, in Hcbrard et al (2002). We performed 
several fits by fixing the column density of a given clement 
X to a different value in each fit. In all of these fits, all 
other parameters are not fixed and we compute the best x 2 
for each value of column density. The plot of A% 2 versus 
N(X) provides the 1, 2, 3, 4, 5cr range for N(X). This is 
illustrated with an example in Figure [HI Examples of the 
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Figure 8. Example of a Ax 2 curve, here for C I in com- 
ponent A. Each point corresponds to a fit performed with 
the C I column density in component A fixed to the value 
given in abscissa. A\ 2 is rescaled by dividing it by the 95 
%-confidence level for this fit (rescaled A\ 2 = 1.41). 3a- 
uncertainty is given by A\ 2 = 3 2 x 1.41 = 12.7. 

fits are presented in Figure Q] for molecules, Figure [2] for 
neutrals, and Figure [3] for ions. 

4.2. Component structure and velocities 

4.2.1. Atoms with ionization potential lower than H i. 

The common asymmetry of the absorption lines of the neu- 
tral species (Figure^]) suggests the existence of two velocity 
components in all of these lines. When fitting simultane- 
ously C i, C i*, C i**, S i and Fe i, a best fit is achieved 
for two components, that we call component A and com- 
ponent B, for the following LSR velocities: Vlsr = 3.8 and 
0.4kms . They may be converted into heliocentric veloc- 
ities by adding 10.2 kms -1 . 

The b-values derived in the fits, &a = 1.5 ± 0.2 and 
6b = 1.6 ± 0.2 kms -1 are not significantly different for the 
three species with different masses. This is strong support 
to the fact that broadening is primarily non-thermal. 

Interestingly, the velocities identified for the two main 
components in the H I Blue Cloud emission spectrum (see 
Table |4]) are identical to those of A and B within the error 
bars. The b- values are also similar to those expected in cool 
gas, with the H I b-values being slightly higher than that 
of S I and C I due to the mass difference. The main com- 
ponents seen in emission are therefore likely to be related 
to the two main components observed in absorption. 

4.2.2. Molecules 

In contrast, the CO lines do not accept a two-component 
solution. The singularity of the CO component is evi- 
denced by the small width of the CO lines, with bco = 
1.6 ±0.2 kms" 1 , similar to that of component A or B alone. 
This width is significantly lower than the velocity differ- 
ence between A and B, 3.3 kms -1 , toward which value bco 
would tend to increase if CO was distributed over the two 
components. When fitting the CO lines in a two-component 
model together with all well-behaved (i.e. not heavily sat- 
urated) neutral lines, unavoidably in the resulting fit most 
of the CO column density is found in one of the compo- 



nents, and the resulting theoretical profile appears slightly 
shifted relative to the observations. This component cannot 
be identified with either A or B, since its velocity shift is 
significantly higher than the uncertainty of the wavelength 
calibration in a single STIS sub-spectrum. Alternatively, 
when fitting the neutrals and CO lines altogether in a one- 
component model, the fit to the neutral lines is significantly 
degraded. The line characteristics (central velocity and line 
width) of CO in emission agree with that derived from UV 
absorption lines. This agreement validates the UV wave- 
length calibration and confirms that CO, unlike H I and 
the atomic species C I and S I, is found at a single velocity. 

We infer that CO is detected primarily in one com- 
ponent at a velocity intermediate to that of A and B, 
VlsrCO = 2.7 kms - . We will call this component C (see 
Figure [T] illustrating the line fit for two CO bands). 

Further information is provided by the CH and CH + 
absorption line profiles. We performed the profile fitting 
of the CH and CH + lines with Owens together with the 
STIS data of the species C i, S i and CO. We made several 
fits with CH and CH + either using one or two components, 
and with a velocity shift as a free parameter to take into ac- 
count uncertainties in the absolute wavelength calibration 
of both the UV and optical spectra. The best results have 
been obtained if (i) CH + is present in both components A 
and B and (ii) most of the CH absorption is in one compo- 
nent, corresponding to the CO component C. This solution 
implies that the optical spectra are shifted by about 0.3 
kms -1 relative to the STIS data, a value consistent with 
the wavelength calibration uncertainty 

We conclude that the bulk of the CO and CH molecules 
are in a single component (C) at a different velocity to the 
bulk of the neutral atoms carbon, sulfur, iron and hydrogen 
and the CH + molecules, which appear to be distributed 
over two well-separated velocity components (A and B). 
Of course it is not excluded that some neutrals are present 
at the velocity of CO, and alternatively that some CO is 
present in the components A and B, at undetectable levels. 
We believe that the three components A, B and C are, in 
fact, physically related. Within this assumption, the three 
velocity components would reflect the particular velocity 
structure associated with local abundance variations within 
a single cloud. 

4.2.3. Ions and neutrals with ionization potential higher 
than H I. 

Almost all ions and neutrals with an ionization potential 
slightly higher than H I: Mg II, Mg I, Mn II, P II, Ni II, 
C II, N I and O I, can be fitted altogether with a unique 
solution for the line of sight structure : component 1 at 2.5± 
0.3 kms -1 , component 2 at — 5.5±0.2kms -1 , component 3 
at -11.0±0.2kms -1 and component 4 at -18.7±0.2kms -1 
(velocities are LSR, see Tabled]). 

The low- velocity component, component 1, corresponds 
to the gas detected in neutral and molecular lines. In the 
unsaturated ion lines this component can be fitted ei- 
ther with two components at the velocity of A and B, 
or with a single component at the intermediate velocity 
Vlsr = 2.5 ± 0.3 kms -1 without affecting the column 
density results or the quality of the fits (monitored by the 
X 2 )- However when adding stronger lines to the fit, it has 
been impossible to find a good fit to the spectra in com- 
ponent 1 with the column density provided by the unsatu- 



rated lines and the A+B velocity structure. This suggests 
that the structure of this component is more complex and 
may indicate the presence of warmer gas at the velocity of 
component 1. Even with a low column density, in strong 
lines warm gas can contribute a large part of the equivalent 
width due to a larger b-value. For the ions, we therefore 
quote a unique column density, obtained using the unsatu- 
rated lines when they exist, for the low velocity component 
1 alone. 

The components structure is summarized in Table [5] 
4.3. Column densities 

4.3.1. Column densities from absorption lines 

Once the velocity structure is defined, column densities can 
be derived following the method described in section 14.11 
In this respect, not all of the lines are equally useful. 

Fe II is the most favorable ion in our sample for two 
reasons. First, several Fe II lines are available with a large 
range of f- values from 0.0014 to 0.32 (see Table \2§. Strong 
lines allow detection of weak components, and weak lines 
enable non-saturated profiles of strong components to be 
analyzed. Another advantage of Fe n lines is that, due to 
its high mass, Fe is less affected by thermal broadening than 
lighter elements; therefore in warm gas Fe n absorption lines 
are narrower and the components better-resolved as seen on 
Figure [31 The errors in N(Fe n) are therefore due only to 
the combination of noise and uncertainty in the continuum 
placement. 

We achieve our poorest quality analyses for ions for 
which only strong lines are available. This is the case here 
for C II, C II* and Si II, for which all components are satu- 
rated and for which we can derive only lower limits. 

When only one weak line for a given element is detected, 
only the strongest component 1 is detected and measured 
and we get column-density upper limits for all other com- 
ponents. This is the case for both Ni n and P II. 

For some elements, we have both a very weak line, 
which is only detected in the strongest component, allowing 
measurement of the column density, and very strong lines 
that are saturated in all components. Combining weak and 
strong lines gives upper and lower limits for the weak com- 
ponents. This is the case for Mg II A1240A and O I A1356A, 
which are respectively 4.4 x 10 4 and 10 3 times weaker than 
the strong lines Mg n A2800A and O i A1302A. S II is an 
intermediate case where the three S n lines are saturated 
in component 1 and at least one S n line is unsaturated 
in the other weaker components. The uncertainty in their 
column density is due to blending between components and 
uncertainty in the continuum placement. 

As mentioned in Section 2.2, Si n* and O I** are de- 
tected primarily at the velocity of component 4, while 
they are hardly detectable in components 1, 2 and 3. 
Nevertheless O I* is detected in the four components. 

Column densities resulting from absorption line fitting 
are provided in Table [HI 

4.3.2. Comparison of absorption/emission column densities 

In this section, we compare column densities derived from 
emission and absorption spectra. 

The total hydrogen column density can be es- 
timated from the color excess, and the standard 



ratio N(B tot )/E( B - V) = 5.8 x 10 21 cm" 2 mag" 1 
(jBohlin et al. 197a ). This provides a measurement of 
N(H to t) = 9.8 x 10 20 cm" 2 . It can also be derived in ad- 
dition using H2 column densities measured in FUSE ab- 
sorption spectra (Table [T]), and the H 1 column density 
estimated from the Blue Cloud H 1 emission spectrum : 
AT(H tot ) = N(R) + 2N(R 2 ) = 9.6 x 10 20 cm- 2 excluding the 
14.5kms _1 background component (Table|l|). The quanti- 
tative agreement between emission and absorption results 
for column densities, velocities, and line widths, supports 
the idea that the two H 1 low velocity components observed 
in the Blue Cloud H 1 emission spectrum correspond to 
components A and B observed in absorption. 

One can also note that there is a consistency between 
the C + column density derived using ISO-LWSdata, and 
the lower limits measured from the saturated absorption 
lines. 

CO column densities were measured both in emission 
and absorption. In absorption, they can be directly com- 
pared to the H2 column densities derived from FUSE 
spectra (Gry et al., 2002). This comparison enables the 
measurement of a total CO abundance N (CO) / N (H2) = 

1- 5t° : 15 io : 7 . 

It was interesting to find that the CO column densities 
were larger when determined in emission than in absorp- 
tion, by a factor 3 in the J = 1, and by more than a factor 
of 7 for the J = 2 level. Because the spatial resolutions 
of the observations are very different, this could imply the 
existence of small-scale structure of the CO emitting gas 
within the SEST beam. On the other hand, part of the 
CO emission could be background emission from the en- 
velope of the nearby Chamaeleon III cloud, which is ob - 
served over the same velocity range (jMizuno et al. 20011 ). 
Based on measurement of the CO excitation, we consider 
that this second possibility is however unlikely. The col- 
umn density ratio between the J = 2 and 1 levels, derived 
from the emission spectrum provides an excitation temper- 
ature T ex = 2.9 K (population ratio N(J = 2)/N(J = 
1) = 0.24), a value very close to the Cosmic Microwave 
Background temperature. The low excitation implies that 
collision excitation of the CO transitions by H or H2, is 
negligible when compared to radiative excitation, hence the 
local density is much lower than the critical density of the 

2- 1 transition nu = 10 4 cm~ 3 . However, if the gas was as- 
sociated with the Cham III cloud, we expect that radiative 
excitation alone could produce a higher T ex . 

4.4. Temperatures 

As noted in Section 14. l\ information on the temperature 
of the gas in each component follows from the ability to 
fit simultaneously lines from several elements of different 
masses, which allows the nature of the broadening to be 
determined. 

4.4.1. Cool gas 

For components A and B, the similarity of the b- values for 
lines from heavy elements of different atomic masses (see 
Section 14. 2p suggests that their widths are dominated by 
non-thermal motions. To estimate the gas temperature, we 
compared heavy-element line-widths to the widths of the 
associated H 1 line emission. 
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Table 5. Interstellar component characteristics derived from the absorption lines in the UV and optical spectrum. The 
ions absorption in component 1 may be distributed over the velocities of components A and B however the fits have been 
performed with one component. Velocities: The listed uncertainties on Vlsr apply to relative velocities. The absolute 
uncertainty from the precision of the STIS wavelength calibration is about 1.2 kms -1 , absolute velocities have however 
been confirmed to about 0.3 kms -1 by the good agreement between component C and the CO emission in the heterodyne 
data, which have an absolute precision of 0.1 kms -1 . Temperature: When possible, temperatures are derived from the 
combination of 6-values from species of different masses. The precise temperatures for components A and B have been 
derived under the assumption that they can be identified to the two components seen in emission in the Blue Cloud H 1 
spectrum. Column densities (in cm -2 ): Numbers in parentheses are powers of 10. All error bars are 3 a uncertainties, 
except for the intervals given by the arrows, indicating and upper limit derived from a faint line and a lower limit derived 
from a strong saturated line. Si 11, C 11 and C 11* profiles arc saturated for all components, implying that only a lower 
limit could be derived for the whole absorption. 



Although the beam of the H 1 emission observa- 
tions is several orders of magnitude larger than the pen- 
cil beam of UV absorption observations (0.5° vs milli- 
arcseconds), we assume that the samplings of the veloc- 
ity field are similar. By computing spectra produced by 
numerical simulations of mildly compressible turbulence, 
Pety and Falgaronc (2000) showed that one-dimensional 
sampling of a turbulent, homogeneous field is similar to 
that performed by an extended beam. This behavior of the 
linewidths is indeed found in mm-line observations of CO 
in emission and absor ption against extragalactic sources 
(jLiszt and Lucas~9 98"). Once the H 1 linewidths arc in- 
cluded in the analysis, the temperatures, yielded formally 
by the 6-value differences between H 1 and heavy neu- 



trals, are 85 ±65 K for A and 270 ±100 K for B. For 
both components, the non-thermal velocity contribution to 
the linewidths is ~ 1.5 km s _1 . We note that these tem- 
peratures apply to the atomic species, and not necessar- 
ily to CO and CH because these molecules appear to be 
kinematically-dccoupled from the atomic species. 

4.4.2. Warm gas 

Components 2 to 4 belong to the second type of interstel- 
lar absorption lines, and T is estimated using the UV-line 
fitting procedure. Component 2 has a temperature typical 
of warm diffuse gas (7 000-10 000 K), while in Component 
3 and 4 the thermal broadening is clearly dominant, with 



temperatures of 13 000±3 000 K for component 3, and about 
20 000 K for component 4 as derived from weak and mod- 
erate lines. In the latter, a more precise determination of T 
is impossible because gas of different temperatures appear 
to be present at this velocity. 

4.4.3. Hot gas 

Indeed the strong Mg II lines show the presence of addi- 
tional high temperature gas at approximately 200 000 K, 
moving at a velocity close to that of component 4. The 
presence of high temperatures is well-illustrated by the clear 
difference between the blue wings of the Mg II profile and 
the N I profile in Figure [3l While the relatively steep blue 
wing of the N I line can be fitted with a temperature close to 
20 000 K for component 4, the extended blue wing in Mg II 
requires the addition of a hot component at the same ve- 
locity, of a temperature ~ 200 000 K, and a column density 
about 10% of that of the total column density of component 
4. If this extended blue wing was a damping wing due to a 
component with very high column density, this component 
would create a significative absorption in the faint Mg II 
line at 1240 A, at least as strong as that of component 1. 
However, no additional component is detected significantly 
in Mg II A1240A. The blue wing can therefore only be due 
to high-temperature, low-column density gas. 

High temperatures are also detected for the excited 
species. The presence of hot gas is traced by the large width 
of the Si II* lines compared to the width of other lines of 
moderate strength. The FWHM of the Si II* line, shown in 
Figure H is ~ 18 kms \ corresponding to ~ 200 000 K. In 
the C II* line, the blue wing is extended, up to LSR veloc- 
ities of — 50 kms - . This may correspond to the presence 
of very high temperature gas at a velocity close to that of 
component 4, as in the case of Mg II. 

The high-ionization species Si III, Si IV and C iv have 
been detected in the spectra at the velocity of component 4. 
They arc broad and the temperature implied by their width 
is at least 200 000 K. A temperature close to 10 6 K, is in 
fact consistent with the observations (see figure 4 for an 
illustration), although the large width of the lines could be 
due to the presence of several components at a tempera- 
ture of a few 10 5 K. OVI is unfortunately not available in 
absorption in the FUSE spectra, because the stellar contin- 
uum is too weak at these wavelengths. 

4.5. Summary of the line of sight 

The description of the line of sight that emerges from the 
direct analysis of the spectral data is complex. The bulk of 
the gas is at low LSR velocity. Neutral low-potential atoms 
and most of the H I arc distributed over two components 
separated by ~ 3 kms -1 (called A and B) that appear in 
these species to have temperatures of at most a few hundred 
K, and are associated with CH + . The molecules CO and CH 
however, are within a single component of similar b-valuc 
(b ~ 1.5 kms" 1 ) at a velocity intermediate between that of 
A and B. 

Other components exist at high negative velocities with 
respect to the bulk of the gas. Negative- velocity gas is seen 
in the absorption lines of ions and atoms with ionization 
potentials higher than that of H I, as well as in emission 
in H I. The temperature of these components ranges from 



typical temperatures of diffuse atomic gas (~ 10 000 K) to 
somewhat higher temperatures (comparable to 20 000 K), 
and increases with velocity offset from the bulk of the gas, 
in addition to the gas excitation state. The highest velocity 
component, while evident for moderately-ionized species, 
of temperature approximately 20 000 K, is associated with 
highly-excited, highly-ionized hot gas of temperatures up 
to several 10 5 K. 

5. The cool gas and its traces of warm gas 

The bulk of the gas detected along the line of sight, is cool 
diffuse molecular gas at low velocities. We comment on its 
properties within the framework of diffuse clouds studies. 

Diffuse molecular clouds are often modeled as extended, 
homogeneous, low-density structures with properties regu- 
lated by the progressive extinction of the UV field by dust. 
In a companion paper (Nehme et al. 2007), we show that 
most observations of the cold, low-velocity gas, namely the 
H2 abundance and ortho-to-para ratio, the CO abundance, 
the C 1 abundance and excitation, and the 158/^m C 11 line 
emission fit with model computations for a diffuse molecu- 
lar cloud of a gas density of 80 cm" 3 , illuminated by a UV 
radiation field close to the mean Solar Neighborhood value. 
The gas temperature, derived from the balance between 
cooling and heating processes, ranges from 60 to 80 K. 
Combining the gas density and gas temperature we calcu- 
late a thermal pressure, P/k ~ 5000cm _3 K, which is within 
the range of values measured by Jenkins (2002), for clouds 
within the local Bubble. While the physical conditions in- 
ferred from the model probably apply correctly to the bulk 
of the gas, we hereafter discuss two points that show the 
limitations of the model's simplicity. 

The first point is the apparent kinematic separation 
of molecular and atomic species in the cool gas, with the 
distribution of atomic species and CH + over two velocity 
components A and B, and the molecules CO and CH ob- 
served at a single intermediate value of velocity. As noted in 
Section 14.21 we believe that these components are related 
and therefore treat them as a single cloud in the model- 
ing. However, the observed velocity structure of the cloud 
cannot be accounted for by the photon-induced chemistry 
model where most of the C 1, CO and H 2 are predicted to 
occupy the same volume (see figures 3 and 4 in Nehme et 
al. 2007). 

The second point is the existence of traces of warm 
molecular g indicated by a CH + absorption line, and a 
population of J > 2 H 2 rotational levels that exceed model 
values by almost two orders of magnitude. Falgarone et 
al. (2005) determined an average Galactic-fraction of warm 
H2, in cool diffuse gas, expressed as the ratio of the col- 
umn density of H2 molecules in levels J > 2 (or A(H 2 ) 
per magnitude of gas sampled: iV(H|)/j4y = 4 x 10 17 
cm" 2 . The warm H 2 predicted along the line of sight to 
HD 102065, with this average Galactic fraction, would then 
equal iV(H 2 ) = 2.7 x 10 17 cm" 2 , compared to the range 1- 
3 xl0 17 cm" 2 given in Table 1. The ratio of the CH+ to 
the total hydrogen column densities towards HD 102065 
N(CR + )/N H = 1.2 x 10" 7 , is comparable to the values 
calculated in previous observations by Crane, Lambert & 
Scheffer (1995) and Gredel (1997). The formation of CH+ 
in diffuse gas involves the endothermic reaction C + + H 2 — > 
CH+ + H (AE/k = 4640 K), and has to be triggered by 
deposition of supra-thermal energy cither in MHD shocks 



(Flower & Pineau des Forets 1998), or in large velocity 
shears (Joulain et al. 1998). 

The chemical patterns triggered by the local deposition 
of non-thermal energy involve several other endothcrmic 
reactions, including CH destruction, which produces neu- 
tral carbon. According to the models explored in Joulain 
et al. (1998), the column density of C I produced in 
the warm chemistry N(C W ), scales with that of CH + as 
N(C W ) pa 37V(CH+) or N(C W ) » 3 x 10 13 cm- 2 . This cor- 
responds to less than 10% of the total C I column den- 
sity iV(C) = 5.7 x 10 14 cm -2 detected in components A 
and B (see Table [5]). The domain of parameters influenc- 
ing the out-of-cquilibrium chemistry is difficult to explore, 
and this is only an estimate provided by existing models. 
The production of CO is enhanced by the warm chemistry 
patterns, due to enhanced production in the warm gas of 
OH, HCO + and CH^~, all chemical precursors of CO. It is 
possible that the observed heterogeneity of CO abundances 
discussed above, originates in such processes. 



6. HD 102065 gas components within the local 
bubble context 

In this section, we place the HD 102065 observations in 
a broader context by relating them to the interaction be- 
tween the nearby interstellar medium, and the Sco-Cen OB 
association, the Local and the Loop I Bubbles. 

Based on an extensive study of colour excess versus dis- 
tance for a large sample of stars within 294° < 1 < 308° 
and -20° < b < 5°, Corradi et al. (1997) concluded that 
the local density of matter is low until an extended in- 
terstellar dust feature at 150 ± 30 pc from the Sun associ- 
ated with the Chamaeleon, Musca and Southern Coalsack 
dark clouds. The same authors complemented their in- 
vestigation with Nal absorption spectra. Absorption lines 
at negative velocities are detected towards many stars in- 
cluding the stars closest to the Sun (Corradi et al. 2004), 
which is indicative of a low column density (Nh ~ a few 
10 19 cm -2 ), nearby (d < 60 pc) sheet of gas, outflowing 
from the Scorpius-Centaurus stars with mean radial ve- 
locity Vlsr = — 7kms~ 1 . The presence of such a sheet 
was first proposed by Egger and Aschenbach (1995) to ac- 
count for the soft X-ray shadow seen in the ROSAT All Sky 
Survey maps. 

In the northern sky, the Loop I shell extends up to North 
Polar radio spur. In the southern sky, a broad H I filament 
at b = —35° and Vlsr ~ — 20kms _1 marks the outer ex- 
tension of gas swept up by the 10-15 Myr old LCC and 
UCL subgroups (de Geus 1992). 

In the Argentine H I survey (Kalberla et al. 2005 
and Bajaja et al. 2005), H I gas at velocities Vlsr < 
— 20kms _1 , down to — 50kms _1 , is visible over 60° in lon- 
gitude about the Sco-Cen sub-groups (see Figure^. The se- 
lection in radial velocity focusses on gas streaming towards 
the Sun at the center of Loop I. We note also that much of 
the negative velocity gas close to the Upper-Scorpius (US) 
group that still contains O stars must be photo-ionized and 
thus not seen in H I. The H I survey data show that neg- 
ative velocity gas is clumped and scattered over a wide 
range of negative velocities, which cannot be accounted for 
within the simple picture of an expanding shell about a lo- 
calized group of massive stars. We propose instead that we 
are looking at dispersed fragments of the parent cloud of 



the LCC and UCL stars accelerated to a range of velocities 
by supernovae blast waves. Since the stellar groups have 
moved with respect to the Sun (de Zeeuw et al. 1999), the 
explosions have occured at different locations and distances 
from the Sun fragmenting and spreading the stars' parent 
cloud between the Local and Loop I bubbles. This view 
accounts for the stream of negative-velocity gas observed 
close to the Sun, including the local cloud about the Sun, 
detected at negative velocities, streaming from the general 
direction of the Sco-Cen association (Ferlet 1999). Our view 
of the relation between the nearby interstellar medium and 
the Sco-Cen stars, the Local and Loop I bubbles builds 
on the original scenario proposed by Breitschwerdt et al. 
(2000). However, the observations suggest a more complex 
picture where there is not a continuous sheet of swept-up 
gas separating the two bubbles, but instead clouds of cold 
and warm gas spread out between the two bubbles. 
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Figure 9. H I emission at negative velocities —50 < 
Vlsr < — 20kms _1 in direction of the Sco-Cen asso- 
ciation. H I observations from Bajaja et al. (2005) and 
Kalberla et al. (2005). The position of HD 102065 and the 
US (Upper-Scorpius), UCL (Upper Centaurus Lupus) and 
LCC (Lower Centaurus Crux) sub-groups of the OB associ- 
ation arc marked. The black area indicates the region where 
the selected range of negative velocities overlaps with the 
Galactic rotational velocities of distant gas. The electronic 
version is in color. 

We believe that the diffuse molecular cloud in the fore- 
ground of HD 102065, is such a cloud. As discussed by 
Mizuno et al. (2001), there is strong evidence that Dcld 
300.2-16.9 (the Blue Cloud) is not part of the Chamaeleon 
molecular clouds complex. A distance of 70 ± 15 pc from 
the Sun, is inferred from the parallax of the T-Tauri star 
T Cha, which coincides in position and velocity (Covino 



et al. 1997) with the cloud CO emission. This places the 
cloud at a distance of ~ 30 pc from the center of the LCC 
sub-group (de Geus 1992). Its cometary structure, with a 
tail extending in the opposite direction to that of the LCC 
stars, suggests that it could have been shaped by a super- 
nova explosion from one of the stars. With a total mass of 
~ 100 M Q (including the tail, as calculated by Boulanger 
et al. (1998) and corrected assuming a distance of 70 pc), 
the cloud is too massive to have been globally accelerated 
by a distant explosion but its structure suggests that gas 
ablated from the surface by the blast wave is streaming out- 
wards. In the H I data, the low velocity components A and 
B are connected to the cloud head but extend in different 
directions along the cloud tail (Figure [5]). This may imply 
that they represent both sides of the gas flow, that they 
encompass the cloud core, and that the difference in radial 
velocity is the result of a difference in angle between the 
streaming direction, and the line of sight. At the surface of 
the gas flow, the CO and CH molecules are photodissoci- 
ated, explaining why they are not detected in components 
A and B, but instead in a single component C, which would 
correspond to the internal part of the gas flow. To account 
for the 3kms _1 velocity difference between components A 
and B within a reasonable angle difference (~ 45°), the gas 
streaming velocity must be ~ lOkms . The tail length 
~ 2° (2.4 pc for a distance of 70 pc) translates into a time 
elapsed since the supernova explosion of ~ 2 — 3 x 10 5 yr. 
The HD 102065 observations were originally motivated by 
the strong mid-IR emission of the foreground cloud inter- 
preted as an unusually high abundance of small dust grains 
(jBoulanger et al. 19901 ). This distinct characteristic may be 
related to the position of the cloud within the Local Bubble, 
and its plausible interaction with a supernova blast wave. 

The interaction between interstellar clouds and shock 
waves travelling in their surrounding medium, has been in- 
vestigated using numerical simulations. Images presented 
in Nakamura et al. (2006) help to picture the complexity of 
the interaction that ablates matter from the cloud and gen- 
erates a turbulent flow streaming away from the cloud. The 
fact that the C IV and Si iv absorption lines peak at neg- 
ative velocities qualitatively agrees with this picture where 
one expects the low density gas in the turbulent mixing 
layer to be entrained by the flow of hot gas. The discrep- 
ancy between gas temperature and ionization state evident 
for some of the Mg II, Si II* and C II* observed at a tem- 
perature 2 x 10 5 K, implies that dynamical mixing of hot 
and warm gas, on timescales shorter than that necessary to 
reach ionization equilibrium, is taking place. 



7. Conclusion 

High resolution spectroscopic observations are presented 
and used to characterize the physical state and velocity 
structure of the multiphase interstellar medium observed 
towards the nearby star HD 102065. These spectroscopic 
observations provide detailed information on the kinemat- 
ics and physical state of the gas, along the line of sight. 

— Gas is observed over a wide range of velocities, with cold 
molecular gas at low velocities, and warmer and lower 
density gas at negative velocities. Most of the hydrogen 
column density is in cold diffuse molecular gas at low 
LSR velocities. 



— The spectra show three distinct components at nega- 
tive velocities as low as Vlsr = — 20kms~ 1 . The ex- 
citation and temperature of the species detected in the 
negative velocity components, increase with the velocity 
separation from that of the diffuse cloud. H I emission 
maps show that the negative velocity gas extends over 
~ 200 deg 2 . It is fragmented and spreads over velocities 
down to — 50kms _1 . 

— A striking result of the data analysis is the detection of 
gas out of equilibrium in both the low and the negative 
velocity gas. At low velocities, we find that the cold 
molecular gas is mixed with traces of warmer molecular 
gas where the observed CH + must be formed and which 
is also traced by the observed excess population of H2 
in the J > 2 levels. At negative velocities, gas in a low 
ionization state (Mg II, Si 11* and C 11*) is observed 
far out of ionization equilibrium at a temperature of 
2 x 10 5 K. 

We have set the observational results within the sce- 
nario proposed by Maiz-Apellaniz (2001) and Bcrghofcr 
and Breitschwcrdt (2002) that connects the origin of the 
Local Bubble to supernovae explosions from early star 
members of the oldest Sco-Cen groups. 

— The gas at low LSR velocities is observed to be asso- 
ciated with the Blue Cloud (Dcld 300.2-16.9) located 
within the local bubble at a distance of ~ 70pc from the 
Sun. Its cometary shape with a tail pointing away from 
the LCC sub-group of the Sco-Cen association ~ 30pc 
away, strongly suggests that it has been shaped by the 
interaction with a recent (2 — 3 x 10 5 yr) supernova ex- 
plosion from one of the LCC stars. The line of sight 
to HD 102065 goes across the cloud tail. The velocity 
structure of the low velocity gas may be understood as 
the two sides of a flow of gas ablated from the cloud head 
Dcld 300.2-16.9. The large abundance of small dust par- 
ticles inferred from the IRAS colors of the Blue Cloud 
tail, possibly results from the cloud interaction with the 
supernova blast wave. 

— The negative velocity gas must have been accelerated 
by a supernova explosion that occurred in the Sco-Cen 
OB association. Our interpretation of the Blue Cloud 
shape would imply that it is physically connected with 
the lower velocity gas. The column density of Si 11* at 
— 20kms~ 1 might be a signature of an ongoing shock. 

A complete diversity of diffuse interstellar-medium compo- 
nents, is present along the line of sight. Our analysis of both 
low and negative velocity, hot, warm and cold gas compo- 
nents, stresses the complex dynamical interplay between 
the interstellar medium phases, in agreement with the pre- 
dictions of numerical simulations (Audit and Henncbcllc 
2005, de Avillcz and Breitschwerdt 2005). We started this 
work by considering an apparently unusual cloud, but now 
believe that the questions we are now asking, have signifi- 
cance for the general study of diffuse, interstellar matter in 
the Galaxy. 
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